Pitch-angle scattering in magnetostatic turbulence. I. Test-particle simulations and the validity of analytical results

Context. Spacecraft observations have motivated the need for a refined description of the phase-space distribution function. Of particular importance is the pitch-angle diffusion coefficient that occurs in the Fokker-Planck transport equation. Aims. Simulations and analytical test-particle theories are compared to verify the diffusion description of particle transport, which does not allow for non-Markovian behavior. Methods. A Monte-Carlo simulation code was used to trace the trajectories of test particles moving in turbulent magnetic fields. From the ensemble average, the pitch-angle Fokker-Planck coefficient is obtained via the mean square displacement. Results. It is shown that, while excellent agreement with analytical theories can be obtained for slab turbulence, considerable deviations are found for isotropic turbulence. In addition, all Fokker-Planck coefficients tend to zero for high time values.Comment: 8 pages, 10 figures, accepted for publication in Astron. Astrophy

Cloud scale influences on mesoscale precipitation patterns

A two dimensional time dependent finite difference grid cloud model is discussed. The model simulates atmospheric motions, potential temperature, water vapor, cloud liquid, cloud ice, rain and small hail. Lateral boundary conditions are open allowing flow in and out of the model domain. Various amounts of convergence were simulated to test the effects on cloud initiation and development. Soundings were run and results discussed

A novel code for numerical 3-D MHD studies of CME expansion

A recent third-order, essentially non-oscillatory central scheme to advance the equations of single-fluid magnetohydrodynamics (MHD) in time has been implemented into a new numerical code. This code operates on a 3-D Cartesian, non-staggered grid, and is able to handle shock-like gradients without producing spurious oscillations. <br><br> To demonstrate the suitability of our code for the simulation of coronal mass ejections (CMEs) and similar heliospheric transients, we present selected results from test cases and perform studies of the solar wind expansion during phases of minimum solar activity. We can demonstrate convergence of the system into a stable Parker-like steady state for both hydrodynamic and MHD winds. The model is subsequently applied to expansion studies of CME-like plasma bubbles, and their evolution is monitored until a stationary state similar to the initial one is achieved. In spite of the model's (current) simplicity, we can confirm the CME's nearly self-similar evolution close to the Sun, thus highlighting the importance of detailed modelling especially at small heliospheric radii. <br><br> Additionally, alternative methods to implement boundary conditions at the coronal base, as well as strategies to ensure a solenoidal magnetic field, are discussed and evaluated

Modulation of the local density of states within the dd-density wave theory in the underdoped cuprates

The low temperature scanning tunneling microscopy spectra in the underdoped regime is analyzed from the perspective of coexisting $d$-density wave and d-wave superconducting states. The calculations are carried out in the presence of a low concentration of unitary impurities and within the framework of the fully self-consistent Bogoliubov-de Gennes theory, which allows local modulations of the magnitude of the order parameters in response to the impurities. Our theory captures the essential aspects of the experiments in the underdoped BSCCO at very low temperatures.Comment: 4 pages, 4 eps figures, RevTex4. New added material as well as reference

Efficient numerical diagonalization of hermitian 3x3 matrices

A very common problem in science is the numerical diagonalization of symmetric or hermitian 3x3 matrices. Since standard "black box" packages may be too inefficient if the number of matrices is large, we study several alternatives. We consider optimized implementations of the Jacobi, QL, and Cuppen algorithms and compare them with an analytical method relying on Cardano's formula for the eigenvalues and on vector cross products for the eigenvectors. Jacobi is the most accurate, but also the slowest method, while QL and Cuppen are good general purpose algorithms. The analytical algorithm outperforms the others by more than a factor of 2, but becomes inaccurate or may even fail completely if the matrix entries differ greatly in magnitude. This can mostly be circumvented by using a hybrid method, which falls back to QL if conditions are such that the analytical calculation might become too inaccurate. For all algorithms, we give an overview of the underlying mathematical ideas, and present detailed benchmark results. C and Fortran implementations of our code are available for download from http://www.mpi-hd.mpg.de/~globes/3x3/ .Comment: 13 pages, no figures, new hybrid algorithm added, matches published version, typo in Eq. (39) corrected; software library available at http://www.mpi-hd.mpg.de/~globes/3x3

Real-Time Operating System/360

RTOS has a cost savings advantage for real-time applications, such as those with random inputs requiring a flexible data routing facility, display systems simplified by a device independent interface language, and complex applications needing added storage protection and data queuing

Sharing emotions and space - empathy as a basis for cooperative spatial interaction

Boukricha H, Nguyen N, Wachsmuth I. Sharing emotions and space - empathy as a basis for cooperative spatial interaction. In: Kopp S, Marsella S, Thorisson K, Vilhjalmsson HH, eds. Proceedings of the 11th International Conference on Intelligent Virtual Agents (IVA 2011). LNAI. Vol 6895. Berlin, Heidelberg: Springer; 2011: 350-362.Empathy is believed to play a major role as a basis for humans’ cooperative behavior. Recent research shows that humans empathize with each other to different degrees depending on several modulation factors including, among others, their social relationships, their mood, and the situational context. In human spatial interaction, partners share and sustain a space that is equally and exclusively reachable to them, the so-called interaction space. In a cooperative interaction scenario of relocating objects in interaction space, we introduce an approach for triggering and modulating a virtual humans cooperative spatial behavior by its degree of empathy with its interaction partner. That is, spatial distances like object distances as well as distances of arm and body movements while relocating objects in interaction space are modulated by the virtual human’s degree of empathy. In this scenario, the virtual human’s empathic emotion is generated as a hypothesis about the partner’s emotional state as related to the physical effort needed to perform a goal directed spatial behavior

Two experiments for the price of one? -- The role of the second oscillation maximum in long baseline neutrino experiments

We investigate the quantitative impact that data from the second oscillation maximum has on the performance of wide band beam neutrino oscillation experiments. We present results for the physics sensitivities to standard three flavor oscillation, as well as results for the sensitivity to non-standard interactions. The quantitative study is performed using an experimental setup similar to the Fermilab to DUSEL Long Baseline Neutrino Experiment (LBNE). We find that, with the single exception of sensitivity to the mass hierarchy, the second maximum plays only a marginal role due to the experimental difficulties to obtain a statistically significant and sufficiently background-free event sample at low energies. This conclusion is valid for both water Cherenkov and liquid argon detectors. Moreover, we confirm that non-standard neutrino interactions are very hard to distinguish experimentally from standard three-flavor effects and can lead to a considerable loss of sensitivity to \theta_{13}, the mass hierarchy and CP violation.Comment: RevTex 4.1, 23 pages, 10 figures; v2: Typos corrected, very minor clarifications; matches published version; v3: Fixed a typo in the first equation in sec. III

Intercomparison of numerical models of flaring coronal loops

The proposed Benchmark Problem consists of an infinitesimal magnetic flux tube containing a low-beta plasma. The field strength is assumed to be so large that the plasma can move only along the flux tube, whose shape remains invariant with time (i.e., the fluid motion is essentially one-dimensional). The flux tube cross section is taken to be constant over its entire length. In planar view the flux tube has a semi-circular shape, symmetric about its midpoint s = s sub max and intersecting the chromosphere-corona interface (CCI) perpendicularly at each foot point. The arc length from the loop apex to the CCI is 10,000 km. The flux tube extends an additional 2000 km below the CCI to include the chromosphere, which initially has a uniform temperature of 8000 K. The temperature at the top of the loop was fixed initially at 2 X 1 million K. The plasma is assumed to be a perfect gas (gamma = 5/3), consisting of pure hydrogen which is considered to be fully ionized at all temperatures. For simplicity, moreover, the electron and ion temperatures are taken to be everywhere equal at all times (corresponding to an artificially enhanced electron-ion collisional coupling). While there was more-or-less unanimous agreement as to certain global properties of the system behavior (peak temperature reached, thermal-wave time scales, etc.), no two groups could claim satisfactory accord when a more detailed comparison of solutions was attempted